2011 Fermi Symposium, Roma., May. 9-12 



1 



o 

(N 

> 
O 

in 



X 

6 



Modeling High-energy and Very-high-energy 7-rays from the Terzan 5 
Cluster 

C. Venter and O.C. de Jager 

Centre for Space Research, North-West University, Potchefstroom Campus, Private Bag X6001, Potchefstroom 
2520, South Africa 
Andreas Kopp 

Institut fur Experimentelle und Angewandte Physik, Christian-Albrechts-Universitatzu Kiel, Leibnizstrasse 1 1, 
241 18 Kiel, Germany 
I. Busching 

Institut fur Theoretische Physik, Lehrstuhl IV: Weltraum- und Astrophysik, Ruhr-Universitat Bochum, 
44780 Bochum, Germany 
A.-C. Clapson 

Max-Planck-lnstitut fur Kernphysik, PO Box 103980, 69029 Heidelberg, Germany 



The Fermi Large Area Telescope (LAT) has recently detected a population of globular clusters (GCs) in high- 
energy (HE) 7-rays. Their spectral properties and energetics are consistent with cumulative emission from a 
population of millisecond pulsars (MSPs) hosted by these clusters. For example, the HE spectra exhibit fairly 
hard power-law indices and cutoffs around a few GeV, typical of pulsed spectra measured for the 7-ray pulsar 
population. The energetics may be used to constrain the number of visible MSPs in the cluster (JV v j s ), assuming 
canonical values for the average 7-ray efficiency and spin-down power. This interpretation is indeed strengthened 
by the fact that the first 7-ray MSP has now been identified in the GC NGC 6624, and this MSP is responsible 
for almost all of the HE emission from this cluster I15I1 . On the other hand, it has been argued that the MSPs 
are also sources of relativistic leptons which may be reaccelerated in shocks originating in collisions of stellar 
winds in the cluster core, and may upscatter bright starlight and cosmic microwave background photons to very 
high energies. Therefore, this unpulsed component may give an independent constraint on the total number 
of MSPs (JVtot) hosted in the GC, for a given cluster magnetic field B and diffusion coefficient ko- Lastly, the 
transport properties of the energetic leptons may be further constrained using multiwavelength data, e.g., to 
infer the radial dependence of fcp and B. We present results on our modeling of the pulsed and unpulsed 7-ray 
fluxes from the GC Terzan 5. 
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1. INTRODUCTION 



The recent Fermi LAT detection of several globular 
clusters (GCs) in high-energy (HE) gamma rays 0, 
B UH, [l6j], very plausibly including Terzan 5, un- 
derlined the importance of modeling collective 7-ray 
emission of millisecond pulsars (MSPs) in GCs. The 
HE spectra are thought to represent the cumulative 
contribution of magnetospheric radiation from a pop- 
ulation of MSPs hosted by the GC. This spectral 
component has been calculated for 47 Tucanae 
and Terzan 5 20] in the context of curvature radia- 
tion (CR) by primary electrons being constrained to 
move along curved magnetic field lines in the magne- 
tospheres of an ensemble of MSPs. 

An alternative calculation [7| considered a scenario 
where inverse Compton (IC) scattering, and not CR, 
was responsible for the HE fluxes seen by Fermi. The 
Fermi fluxes may be reproduced for certain parame- 
ters, and this model also predicted very-high-energy 
(VHE) components in some cases. The discovery 
of the luminous MSP PSR J1823-3021A in the GC 
NGC 6624 [l5[ however implies that such putative 
unpulsed HE IC components may be dominated by 
a pulsed CR component, at least for this particular 
GC. 

In addition to the pulsed flux, a steady flux from 



GCs is also expected in the VHE domain. The model 
of @ predicted HE and VHE fluxes from GC by con- 
sidering relativistic leptons escaping from the embed- 
ded MSP population and upscattering soft photons 
from background radiation fields via the IC process. 
Indeed, some MSPs may produce leptons with TeV 
energies due to acceleration of these particles by very 
large magnetospheric electric fields |6| . These leptons 
may be further accelerated in shocks in the GC result- 
ing from colliding pulsar winds. 

A similar calculation of the unpulsed IC component 
to that of Q was performed [TH, [2(| using a parti- 
cle injection spectrum calculated from first principles 
and which is the result of acceleration and CR losses 
occurring in the MSP magnetospheres. No further 
particle acceleration was assumed after escape from 
the magnetosphere, yielding predictions that should 
be considered as lower limits for the VHE flux band. 
This model predicted that 47 Tucanae and Terzan 5 
may be visible for H.E.S.S., depending on the assumed 
model parameters, particularly TVtot and cluster mag- 
netic field B. This model furthermore fixed the par- 
ticle efficiency r/ e to ~ 7% of the average spin-down 
luminosity, reducing the number of free parameters. 

Recent H.E.S.S. upper limits on the TeV emission 
from the 47 Tucanae ji] implied that AT tot - 30 - 40 
for B ~ 10/iG, but N tot becoming quite larger for 
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Figure 1: Differential CR and IC spectra for Terzan 5. 
The green lines are the CR spectrum and errors 
from [20|], scaled using iV v i s = 60 to fit the Fermi data; 
the magenta lines are IC models from Q; Fermi LAT 
data [2| are plotted in black; the black line is an IC 
prediction from [5(, scaled to the new distance and 
bolemetric luminosity, average spin-down luminosity of 
1.8 x 10 34 erg s _1 , and iVtot = 60; the blue lines represent 
IC spectra for different values of B (this work); the red 
line is the H.E.S.S. sensitivity. 



B < 5^G or B > 30^G [H. Also, the Fermi LAT 
HE spectrum implied that there are iV v i s ~ 50 — 60 
MSPs in the cluster (0] inferred N vis - 50). 



2. MODEL 

We have previously calculated the pulsed CR spec- 
trum resulting from relativistic leptons which are ac- 
celerated in the MSP magnetospheres by large pair- 
starved electric fields [l9[ as predicted by the pair- 
starved polar cap (PSPC) model [lH, prior to the de- 
tection of the HE spectrum by Fermi [2] ■ This calcu- 
lation may now be scaled to the Fermi data to infer 
7V vis (Section©. 

We also calculate an unpulsed IC flux component, 
using a cumulative injection spectrum made up of 
electrons leaving the MSP magnetospheres after hav- 
ing been accelerated by these pair-starved magneto- 
spheric electric fields, and neglecting reacceleration in 
the GC. Using updated structural parameters, dis- 
tance, and a much larger bolometric luminosity, we 
calculate the radiation losses and resulting unpulsed 
IC fluxes shown in Fig. [TJ assuming Bohm diffusion 
and GC magnetic fields of B = 1 fiG and B = 10/iG 
(for more details, see pH, We used two radi- 

ation zones in the cluster: a core region extending 
from r — up to r c , as well as a halo region extending 
from rc up to r t . The steady-state particle spectrum 
was approximated by the product of the particle in- 
jection spectrum, and an effective times cale taking 
into account the IC, SR and particle escape loss time 
scales. 

We used both bright starlight and CMB as soft pho- 
ton targets in our IC calculation. The energy density 
of the first was assumed to be ~ 2.4 x 10 4 eV cm -3 
and ~ 1.6 x 10 3 eV cm~ 3 for each of the regions (for 
a temperature T = 4, 500 K, and due to large stellar 
luminosity and small core radius), while the energy 
density for the CMB was taken to be ~ 0.27 eV cm~ 3 
(for T = 2.76 K). 



H.E.S.S. recently detected a VHE excess in the di- 
rection of Terzan 5 offset from the center of the 
GC by 4', and having a size of 9.6' x 1.8' (compared to 
the Fermi maximum likelihood source position which 
is offset from the GC center by 2.4', still within the 
95% source position uncertainty of r g5 = 2.9', and 
source^ extent of 9'). In addition, diffuse X-ray emis- 
sion Q, as well as several radio structures [|| have 
been measured from this GC. There have also been 
updated measurements of Terzan 5's distance (d = 
5.9 ± 0.5 kpc) [H H3, core radius (r c = 0.15'), half- 
mass radius (fhm = 0.52'), tidal radius (r t = 4.6'), 
and total luminosity (L~8x 10 5 L Q ) [l3j . 

In this paper, we present pulsed CR as well as un- 
pulsed IC calculations for an ensemble of MSPs in the 
GC Terzan 5. Independent constraints may be derived 
on iVtot using both of these components, while B may 
be constrained using synchrotron radiation (SR) and 
IC flux components. 



3. RESULTS 

Figure [1] shows the differential CR and IC spectra 
calculated for Terzan 5. Scaling the pulsed CR com- 
ponent [2(| to fit the Fermi LAT data [|J implies that 
the number of visible MSPs is about iV v i s ~ 60 ± 30. 
This finding is consistent with the estimate of 180^io 
obtained by Q , and formally presents a lower limit to 
iVtot- However, the pair-starved model probably over- 
predicts the CR flux by a factor of a few, and further- 
more may not be valid for all MSPs (as inferred from 
light curve modeling (2l| which imply copious pair 
creation in some of the MSP magnetospheres, despite 
their low magnetic dipole surface fields), so that iV v i s 
may be even larger than ~ 60. The CR spectrum 
furthermore cuts off below the H.E.S.S. sensitivity, so 
that the VHE signal probably originates from IC scat- 
tering. 

The VHE fluxes in Figure [1] include IC models 
from [7], scaled predictions from as well as our 
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Figure 2: Constraints on iVtot vs - GC cluster magnetic 
field B. Green line: ICS component using H.E.S.S. 
sensitivity (this work; light green area is excluded). Red 
and blue lines: 7V v i s from 0] and this work (CR 
component); magenta and black lines: number of radio 
MSPs detected and estimated [Tl| . Errors not shown. 



ing this model, we could constrain TVtot and B. Our 
unpulsed spectral results should be re-assessed in view 
of the recent H.E.S.S. detection of Terzan 5 in order 
to obtain updated constraints on these parameters. In 
particular, the observed spectral shape implies that 
reacceleration of particles may be taking place within 
the GC. The offset nature of the source with respect 
to the GC center provides a further puzzle, because if 
the MSPs are located within the GC core radius, the 
HE and VHE emission should be GC-centered. Lastly, 
availability of multiwavelength data on this GC may 
allow us to constrain the radial profile of the soft pho- 
ton energy density, diffusion coefficient ko and cluster 
field B in future. 
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calculated IC spectra for different values of B. While 
the prediction corresponding to B = 1 /xG is just be- 
low the H.E.S.S. sensitivity, the one corresponding to 
B = 10 /iG is above this detection threshold above 
several TeV (assuming N tot = Avis = 60). 

Figure [2] indicates constraints on iVtot vs. B. The 
green line indicates the constraint derived from our 
ICS calculation using H.E.S.S. sensitivity (the light 
green area is excluded). The red and blue lines are 
constraints on iV v i s from Fermi [2[ and our CR calcu- 
lation, while the magenta and black lines indicate the 
number of detected and estimated radio MSPs 



11] 



The halo size was terminated at r t , assuming that 
the soft photon energy density is sufficiently low out- 
side of the GC that IC emission beyond rt may be 
neglected. However, in view of the new HE and VHE 
data that imply that this source is extended in gamma 
rays, one may need to reassess this assumption. In- 
deed, an alternative model found [7| that most HE 
emission may come from a region outside of the GC 
core beyond a radius of 10 pc. 



4. CONCLUSION 

We have obtained pulsed and unpulsed fluxes from 
the GC Terzan 5, assuming CR and IC processes in- 
volving TeV leptons from a number of host MSPs. Us- 
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